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1. Introduction 

Over the last decade, transport properties of one-dimensional systems consisting of 
coupled oscillators on a lattice have been the subject of many theoretical and numerical 
studies, see the review papers [HI [HI EH]- Despite many efforts, our knowledge of the 
fundamental mechanisms necessary and/or sufficient to have a normal diffusion remains 
very limited. Nevertheless, it has been recognized that conservation of momentum plays 
a major role and numerical simulations provide a strong evidence of the fact that one 
dimensional chains of anharmonic oscillators conserving momentum are super-diffusive. 

In this paper we propose a new class of models for which anomalous diffusion 
is observed. The system under investigation presents several analogies with standard 
chains of oscillators, but, and it is our main motivation, has a simpler mathematical 
structure. 

1.1. Description of the model 

Let U and V be two potentials on M and consider the Hamiltonian system {uj(t))t>o = 
( r(t),p{t) )t>Q described by the equations of motion 

^ = V\r,+,)-V\r,), ±^ = U\p,)-U'{p,+^), x G Z, (1) 
at at 

where px is the momentum of particle x, qx its position and = qx — qx-i the 

"deformation". Standard chains of oscillators are recovered for a quadratic kinetic 

energy U{p) = p'^/2. The dynamics conserves three physical quantities: the total 

momentum ^^.Pa;, the total deformation ^^.r^,. and the total energy J2x^x with 

Sx = Vi^x) + U{px). Consequently, every product probability measure z^^,a,a' defined by 

dvux^'iv) = Yl ^">~^ {-(^^^ - - ^''^x} dvx dpx (2) 

is invariant under the evolution. For later purposes, let us also introduce the short-hand 
notation 
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In this paper, we are mainly interested in the case U = V (except for instance in 
Section |44I) . which has the advantage of introducing more symmetries into the problem. 
The state of the system at time t is then more conveniently described by the variable 
rjit) = {rix(t);x G Z} G with ?72x-i = i^x and ri2x = Px- The dynamics can be 
rewritten as: 

dr]x{t) = (V'{7]x+i) - V'{r]x--,)yt. (3) 

The system can therefore also be interpreted as a fluctuating interface where the 
algebraic volume of the interface at site x G Z is given by rjx and the energy by V{rix). 
We focus our study on the (anomalous) diffusion of the energy. The three quantities, 
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momentum rj2x, deformation ^^g^ ?72x+i and energy ^^^^ V{rjx), are conserved but, 
for reasons which will become clear later, our interest lies only in the "volume" ^^.g^ 1x 
and in the energy. 

1.2. Hydrodynamic limit 

Energy transport properties depend strongly on the chosen time-scale. The first natural 
scale to consider is the hyperbolic scale where the system is followed on long times Nt, 
the space being renormalized by a factor A^~^, with N — )■ oo. 

To this end, we define the empirical energy/momentum/deformation measure as 
follows: 

7r^(t, dq) = Yl ^-(^^) ^-/N{dq). q G M- (4) 

At time t = 0, this measure is supposed to converge in probability to some macroscopic 
profile Uolq) dq which has a density w.r.t. the Lebesgue measure. If we assume 
that a local equilibrium hypothesis holds, it is not difficult to show that the expected 
macroscopic evolution equation for Il(t,q)dq = lirajy^oo '^'^ {t,dq) is given by a triplet 
of compressible Euler equations: 

dtii + d,i{u) = o, n(o,-) = no(-), (5) 

where /(H) G M.^ is the macroscopic current whose explicit expression is not important 
here. Unfortunately, proving the local equilibrium hypothesis for ([T]) is out of the 
range of the current mathematical techniques. The main difficulty comes from our 
inability to show that the dynamics is ergodic (in the sense of Definition [1] below). 
If the ergodicity is proved, in the time interval where ([5]) has a smooth solution, the 
relative entropy method of Yau (see |13]) can be adapted to show that the system 
has ([5]) as a hydrodynamic limit. Hence, the problem can be reduced to proving the 
ergodic behavior of the underlying dynamics. Deriving the convergence to ([5]) after the 
shocks is considerably more difficult since even the concept of a solution to a system of 
conservation laws is not fully understood. 

To overcome the above mentioned lack of ergodicity of deterministic systems, it 
has been proposed to add a stochastic perturbation to the dynamics. The theory 
of stochastic perturbations of Hamiltonian dynamics has a long history. To our 
knowledge, the first paper on the ergodicity of infinite lattice models is [H] (see also 
P, El [16], [m [m [26], [32]). The added noise must be carefully chosen in order not to 
destroy the conservation laws we are interested in. In the general case U V, the 
Hamiltonian dynamics can be perturbed by a local noise acting on the velocities (as 
proposed in [16]) but conserving the three physical invariants mentioned at the end of 
Section [LT] (see Section [131). With such additional noises, the perturbed dynamics can 
be proved to be ergodic (see Theorem [3l), so that ([5]) is obtained as a hydrodynamic 
limit. However, our motivation being to simplify as much as possible the dynamics 
considered in [H [2] without destroyin g the anomalous behavior of the energy diffusion. 
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we mainly focus on the symmetric case U = V with a noise conserving only the two 
important quantities responsible of the anomalous transport behavior, namely the the 
energy and the volume. Thus, we introduce a new stochastic energy-volume conserving 
dynamics, which is still described by ([3]) between random exponential times where two 
nearest neighbors heights ijx and rj^+i are exchanged (see Subsection 12.31 for a precise 
definition). Observe that the noise still conserves the total energy and the total volume 
but destroys the conservation of momentum and deformation. Therefore, only two 
quantities are conserved and the invariant Gibbs measures of the stochastic dynamics 
correspond to the choice A = A' in ([2]). We denote z^/3,a,a (resp. Z{f3, A, A)) by /i/j^A (resp. 
Z{P, A)) and we use in the sequel the short-hand notation 

The first main result of this paper is that the perturbed dynamics informally 
described above is ergodic (see Theorem [2]). Consequently, before the appearance of 
shocks, the stochastic energy-volume conserving dynamics has a hyperbolic system of 
two conservation laws as a hydrodynamic limit (see Theorem [1]). 

1.3. (Super) Diffusive limit 

In the second part of the paper we investigate the diffusion of the energy at a longer time- 
scale. If the process has a diffusive behavior then the relevant time scale is the diffusive 
one, where the system is studied over long times N'^t with space renormalized by a 
factor N~^. We claim that the system genuinely displays an anomalous energy diffusion 
so that the diffusive scale is not the relevant one. Heuristically, we can interpret this 
anomalous diffusion as a consequence of the volume conservation law (see Section [2^2]) . 

We start the infinite system ([T]) under the equilibrium distribution fi/s^y 
consider first the fluctuation field in the hyperbolic time scaling: 

yNif G) = ^ E G (^) ® ie.m - e) , (7) 



where 

( GM \ 
Giy) = GM 
\ GM ) 

is a smooth vector valued test function with compact support. We expect that 3^Ar(t, ■) 
converges in law to y{t, ■), where y is solution of the linearized equation 

dty + DI{e)d,y = Q, (8) 

with DI{e) the differential of I at e. Hence, in the hyperbolic scaling, fluctuations 
evolve deterministically according to 

3^(t,G) = 3^(0,e*"*G), 
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where U = [DI{e)]dg and U* = —[DI{e)]*dg with [DI{e)]* the transpose matrix of 
DI{e). 



To see a nontrivial behavior of the fluctuation field, we need to look at on a 
longer time scale tN^~^°', for some a > 0. It is expected that, after subtracting the 
transport term appearing in the hyperbolic time scale, the field 



converges to some limiting field y. The case a = 1 would correspond to a diffusive 
behavior with 3^ the solution of the linear stochastic partial differential equation 



where W{x,t) is a standard space-time white noise. Here, x is the compressibility and 
X>°° = limt_>.oo T^i3,\,y {t) the limiting dijfusivity (see (jHJ and fH5]l for the definitions of 
these quantities). We refer the reader to [39] for a general background reference, and to 
[27] for a rigorous proof of the convergence for asymmetric simple exclusion processes. 

For the class of models we consider, our conjecture is that a is in general strictly 
lower than 1. The value of a and the nature of y are not expected to be universal and 
should depend on some specific properties of the potentials. We also expect a similar 
picture when the deterministic dynamics ([1]) is replaced by the stochastic energy- volume 
conserving dynamics (we then denote by I^/3,a and Cjs^x the corresponding diffusivity and 
current-current correlation function). This anomalous value of a should be reflected in 
the divergence of the diffusivity Vj^^x^yit) {Vj^^xit) for the stochastic dynamics) in the 
large time limit t — ?■ oo. 

We are not able to study theoretically this problem for the deterministic dynamics 
and we have to turn to computer simulations of nonequilibrium systems in their steady- 
states. A chain of length 2A^+1 is coupled at each extremity (left and right) to a thermal 
reservoir fixing the temperature (T^ on the left, on the right). In the stationary state 
(■)ss, the average current {Jn)ss is measured (see Section [6] for more precise definitions). 
The quantity of interest is the divergence exponent 5 of the transport coefficient 



Anomalous diffusion corresponds to 5 > 0. 

For a normal transport, the link between the two situations can be seen through a 
Green-Kubo formula for the limiting diffusivity, which expresses the latter as a quantity 
proportional to the time integral of the equilibrium current-current correlation function 
(see (Ho!) and (148!) for a precise definition). It is widely accepted, but not proved, that for 
normal diffusive systems the limiting diffusivity coincides with the transport coefficient 
K = YmvN^oo i^N ■ For anomalous diffusion transport, k = +oo and the current-current 
autocorrelation function which appears in the Green-Kubo formula is not integrable 
because it decays too slowly. Consequently the limiting diffusivity is infinite. 

Our second main results are the following. First, we show numerically that, for 
generic anharmonic potentials V , the dynamics (JSj) has an anomalous diffusion. We 




(9) 
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also show that this phenomenon persists if the deterministic dynamics is replaced by 
the stochastic energy-volume conserving dynamics, and that the divergence exponent S 
depends on the strength of the random perturbation. Secondly, for the stochastic 
energy- volume conserving dynamics and a harmonic potential V{r) = we compute 
explicitly the equilibrium current- current correlation function Cp^\{t) and show that it 
decays as for large t (see Theorem H]). This implies that the diffusivity T)p x{t) 
diverges as which is a clear manifestation of the super- diffusion of the energy for 
this model. 

1.4- Organization of the paper 

We present more precisely the model under investigation in Section [21 We first study the 
hydrodynamic limit in Sections |3] (presentation of the general result) and H] (proof of the 
fundamental ingredient for the limit to hold, namely the ergodicity of the dynamics). We 
then consider diffusion properties, starting with analytical results on the longtime tail 
of the current autocorrelation function, which can be obtained for harmonic potentials 
(see Section [5]), and providing then scalings of the energy current obtained by numerical 
simulations of nonequilibrium systems in their steady-states in Section [61 Some proofs 
are gathered in the Appendix. 

2. The models 

2.1. The deterministic models 

2.1.1. Finite systems Consider the finite box A^r = {— A^, . . . , A^} C Z (with A^ > 1). 
The product space M^^ is denoted by ^at, and a typical element of f^Ar is 77 = {rjx G 
M; X G Ajv}. The deterministic finite volume dynamics {rj^ {t))t>Q G Vt^ is defined by 
its generator 

N-l 

An= {^'(v.+i) - V'{r].^,))d,^ - V'ivN-i) d,^ + \^'(r/-,v+i) d,_,, (10) 

where \^ is a smooth convex potential such that the partition function 

/oo 
exp i-l3V{r) - Xr) dr 
-00 

is well defined for /3 > and A G M. The following microscopic energy-volume 
conservation laws hold for x = —N + 1, . . . , A^ — 1: 

where V is the discrete gradient defined, for any function m : Z — M, by (Vm)(x) = 
u{x + 1) —u{x), and where the microscopic energy and volume currents are respectively 

fx,x+i = -V'ivx+i)V'ivx), Jlx^i = -{V'ivx) + V\vx+^)). 
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2.1.2. Infinite systems. The dynamics in the infinite volume A = Z, with formal 
generator 



is also very important. Since the state space is unbounded, explosion problems can 
arise and the construction of solutions of the dynamics in infinite volume may become 
a technically non trivial problem. To avoid such issues, we restrict ourselves (apart 
from Section [6]) to the case <V" < C for some positive constant C > 0. Then, the 
construction is quite standard since the pioneering work [25] of Lanford et al. We refer 
the interested reader to [IS [16] for further precisions. For any a > 0, let be the set 
of configurations rj such that 

and equip Vl = na>o^a with its natural product topology and its Borel cr-field. The set 
of Borel probability measures on Q is denoted by V{Q). A function / : f2 — i- M is said to 
be local if it depends of t] only through the coordinates {rix ; x G Aj}, Af being a finite 
box of Tj. We also introduce the sets Cq{Q) {k > 1) of bounded local functions on Q 
which are differentiable up to order k with bounded partial derivatives. 

For each initial condition a G f2 the existence and uniqueness of a solution to ([3]) can 
be proved by a classical fixed-point argument a la Picard. The solution ?7(-) := ri{-,a) 
defines a process with continuous trajectories. Moreover each path ri{-, a) is a continuous 
and differentiable function of the initial data a. We define the corresponding semigroup 
{Pt)t>o by {Ptf){cr) = f{rj{t,a)) for any bounded measurable function / on f2. The 
differentiability with respect to initial conditions shows that the Chapman- Kolmogorov 
equations 



(P,/)(a) = /(a) + f\APsf){cr) ds, f G C',{n) 
Jo 

{PJ){a) = f{a) + / {PsAf){a) ds, f G C'o{n). 
Jo 



are valid. With these equations, probability measures u G V{Q) invariant by the 
deterministic dynamics are characterized by the stationary Kolmogorov equation: 



V/GCo^(fi), J {Af){v)du{r^) = 0. 



Denoting the usual scalar product between two vectors a, 6 G by a ■ 6, and recalling 
the notation introduced in ([6]), it is easily seen that every product measure fijs^x defined 
by 

dfifi,x{r]) = n ^^^^ "L-^ ■ ^x} dr]^ 

is invariant for the infinite dynamics. In the sequel, we denote the average of a function 
/ with respect to fif^^x by (/)/3,a- 
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Depending on the potential V at hand, the properties of the dynamics can be 
very different. In the next subsection we discuss the case V{r) = which leads 

to a linear dynamics. A second remarkable potential is the exponential potential 
Vkvm{(i) = e"'' + g — 1, corresponding to the so-called Kac-van-Moerbecke system, 
which is integrable (see [23]). The corresponding system is related to the famous 
Toda lattice [H], i.e. a chain of oscillators with coupling potential Vkym, by a simple 
transformation. 



2.2. The deterministic linear model 

We consider here the specific case V{r) = The dynamics is then linear and can 

be solved analytically using Fourier transform. To simplify the exposition we consider 
the dynamics in infinite volume. We introduce the kth mode fi[k, ■) for k E T = M/Z, 
the one-dimensional torus of length 1 : 

rj{t,k) = Y,V,{t) e^^-^^ 

Then, the equations of motion are equivalent to the following decoupled system of first 
order differential equations: 

^(t, k)=iuj{k)fi{t, k), 
dt 

where the dispersion relation u{k) reads 

u{k) = — 2 sin(27rA;), 
and the group velocity Vg is 

Vg{k) = uj'{k) = —4:71 cos {271 k). 
By inverting the Fourier transform, the solution can be written as 

r/,.(t)= [ riit,k)e-^'^''^dk. 



Note also that the energy of the kth mode 

Ek{t) = ^m,k)\^ = Ek{0) 

is conserved by the time evolution, and that the total energy current = Xl^ez-^ 
takes the simple form 



e 

x,x+l 



r = I Vg{k)Ekdk. 
Jt 



We interpret the waves r/(/c, t) as fictitious particles similar to phonons in solid state 
physics. In the absence of nonlinearities, they travel the chain without scattering. If 
the potential is non-quadratic, it may be expected that the nonlinearities produce a 
scattering responsible for the diffusion of the energy. Nevertheless, the conservation of 
the volume, which is expressed by 

r/(t,0) =7/(0,0), (12) 



9 



plays a crucial role. The identity fll2p is valid even if V is not quadratic. It means 
that the 0th mode is not scattered at all and crosses the chain ballistically. In fact, 
the modes with small wave number k do not experience a strong scattering and they 
therefore contribute to the observed anomalous diffusion of energy. For anharmonic 
chains of oscillators, a similar picture arises with r/(/c) replaced by the phonons. As 
the conservation of momentum for these chains is responsible for the small scattering 
of phonons with small wave numbers, the conservation of the volume for the model 
considered in this paper is responsible for the small scattering of the waves r]{k) with 
small wave numbers. 

2.3. Stochastic energy-volume conserving dynamics 

We now consider energy- volume conserving stochastic perturbations of the deterministic 
dynamics generated by ^ or ^jy- The generator of the finite volume perturbed dynamics 
is written as 

Cn = An + iSn, (13) 
where Sn is the generator of the noise and 7 > its intensity. The generator Sn reads 

N-l 

{SNf){v)= E [fiv^'^^')-fiv)], (14) 

x=-N 

where ?^^'^+^ is the configuration obtained from 77 by exchanging the variables rj^ and 

Vx+l- 

By arguments similar to the one used to prove the well-posedness of the infinite 
deterministic dynamics, it can be shown that the stochastic energy-volume conserving 
dynamics in infinite volume is also well defined (see [16] for details). Its formal generator 
is given by £ = ^ + 7^ where 

(5/)(r/) = $^[/(r/^'^+^)-/(r/)], 

x&Z 

and the corresponding Chapman-Kolmogorov equations (fTT]) are valid for this process 
upon replacing A by C In particular the probability measures /x^^a are still invariant. 

3. Hyperbolic scaling 

We present in this section the hydrodynamic limit of the models described in the previous 
section. To this end, we first need to define some thermodynamic quantities useful 
to describe local equilibria (Section 13. ip . We then informally describe the expected 
hydrodynamic limit in Section 13.21 and conclude this section by stating precisely the 
convergence result (Theorem [1]) in Section [3^ 
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3.1. Definition of thermodynamic variables 

Recall that the probability measures /i^^a form a family of invariant probability measures 
for the infinite dynamics defined in Section O The following thermodynamic relations 
(which are valid since we assumed that the partition function Z is well defined on 
(0, +oo) X M) relate the chemical potentials /3, A to the mean volume v and the mean 
energy e under 

^(/3, A) = /i/3,A(r/.) = -9a ( log Z(/3, A)) , 

/ \ ^^^^ 

e(/3,A) =/i^,A(r(r/,.)) = -dp(\ogZ{l3A))- 

These relations can be inverted by a Legendre transform to express /3 and A as a function 
of e and v. Define the thermodynamic entropy 5" : (0, +oo) x M — t- [— oo, +oo) as 



S{e,v)= inf |/3e + At; + log Z(/3, A)|. 



;,/3>o 

Let U be the convex domain of (0, +oo) x M where S'(e,w) > — oo and U its interior. 
Then, for any (e, := (e(/3, A), w(/3, A)) G lA, the parameters /3, A can be obtained as 

P = {d,S){e,v), X = {d,S){e,v). (16) 

We also introduce the tension r(/3, A) = Ai^,A(^'('7o)) = —^/P- Then, 

^^pAfx,x+l) = ^^P^Ul,x+l) = -2^- (17) 

In the sequel, with a slight abuse of notation, we also write r for r(/3(e, w), A(e, v)) where 
/9(e,f) and A(e,w) are defined by the relations (fT6|) . 

c?.^. Description of the hydrodynamic limit 

Consider the finite closed stochastic energy-volume dynamics with periodic boundary 
conditions, that is the dynamics generated by £Ar,pcr = -^Ar.pcr + iSN^^er where 



(^iV.per/) (t/) = [V'{Vx+l)-V'{n,.{)]dr,J{7l)^ 



and 



with Tn = ]R/(A^Z) is the discrete torus of length A^. We choose to consider the 
dynamics on Tat rather than on Z to avoid nontrivial technicalities. We are interested 
in the macroscopic behavior of the two conserved quantities on a macroscopic time-scale 
Nt as N ^ oo. 

We assume that the system is initially distributed according to a local Gibbs 
equilibrium state corresponding to a given energy-volume profile Xq : T — t- W: 
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X{t, 



in the sense that, for a given system size A^, the initial state of the system is described 
by the following product probability measure: 

where {f3o{x / N) , \o{x / N)) is actually a function of {to{x / N) , \Jo{x / N)) through the 
relations f fT6|) . 

Starting from such a state, we expect the state of the system at time t to be close, 
in a suitable sense, to a local Gibbs equilibrium measure corresponding to an energy- 
volume profile 

' e(t,-; 

satisfying a suitable partial differential equation with initial condition Xq at time t = 0. 
In view of f|T7|) . and assuming local equilibrium, it is not difficult to show that the 
expected partial differential equation is the following system of two conservation laws: 

- oy = 0, 

' (20) 
dt^ - 2dgT = 0, 

with initial conditions e(0, ■) = Cq{-), 0(0, ■) = tio(')- We write ( 120|) more compactly as 
dtX + dgZ{X)=0, X{0,-)=Xoi-), 

with 

The system of conservation laws (120|) has other nontrivial conservation laws. In 
particular, the thermodynamic entropy S is conserved along a smooth solution of (120]) : 

dtS{t,v) = 0. (22) 

Since the thermodynamic entropy is a strictly concave function on U, the system ( l20l) 
is strictly hyperbolic on W (see [37]). The two real eigenvalues of {DZ)(^) are and 
— [i9e(T^) + 29^ (r)], corresponding respectively to the two eigenvectors 

(23) 

It is well known that classical solutions to systems of > 1 conservation laws 
develop shocks, even when starting from smooth initial conditions. Nevertheless, the 
Cauchy problem is locally well-posed in the Sobolev spaces H'^{T) (for s > 3/2). If 
we consider weak solutions rather than classical solutions, then a criterion is needed to 
select a unique, relevant solution among the weak ones. For scalar conservation laws 
(n = 1), this criterion is furnished by the so-called entropy inequality and existence 
and uniqueness of solutions is fully understood. If n > 2, only partial results exist 
(see [371 E])- This motivates the fact that we restrict our analysis to smooth solutions 
before the appearance of shocks. 
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3. 3. Derivation of the hydrodynamic limit 

We now turn to the question of deriving the system of conservation laws as the 
hydrodynamic hmit of the interacting particle system under investigation. We assume 
that the potential V satisfies the following 

Assumption 1. The potential V is a smooth, non-negative function such that the 
partition function Z{f3, A) = exp {—f3V{r) — Ar) dr is well defined for f3 > and 
A G M and there exists a positive constant C such that 



and 



< V"{r) < C, (24) 

lim sup G (0, +oo) , (25) 

\V'(r)V 

lim sup ^ , < +00. (26) 

The hypothesis f l24|) allows to define easily the dynamics in infinite volume; (l25l) 
is needed in the proof of Theorem [21 ( 126|) ensures that the currents of the conserved 
quantities are bounded by the energy. This is useful to introduce a suitable cutoff for 
the derivation of hydrodynamic limits (see [U Section 3]). 

Provided we can prove that the infinite volume dynamics is ergodic in a suitable 
sense (see Definition [1] below), then we can rigorously prove, using the relative entropy 
method of Yau, that ([5]) is indeed the hydrodynamic limit in the smooth regime, i.e. for 
times t up to the appearance of the first shock (see for example [211 112]). Recall indeed 
that a simple computation shows that the stochastic perturbation does not modify the 
hydrodynamic limit since the effect of the latter is observed in the diffusive scale only. 

In most cases, the derivation of hydrodynamic limits is performed for stochastic 
interacting particle systems which are trivially ergodic by construction. For 
deterministic systems on the other hand, the ergodicity is extremely difficult to prove. 
We are only able to show a weaker form of such ergodicity for the process generated by 
A (with the additional assumption that the invariant measure is exchangeable). This 
weaker form is nonetheless sufficient to show that the process generated by C is ergodic 
(see Theorem [2] below). 

As argued in [12], it turns out that when there are more than one conservation 
laws, the conservation of thermodynamic entropy (l22il must hold for the hydrodynamic 
limit to be well defined. This relation is indeed fundamental for Yau's method where, 
in the expansion of the time derivative of relative entropy, the cancelation of the linear 
terms is a consequence of the preservation of the thermodynamic entropy. 

Averages with respect to the empirical energy-volume measure are defined, for 
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continuous functions G, if : T — )■ M, as 
SN{t,G) ' 



VN{t,H) 



1 / X 



We can then state the following result. 

Theorem 1. Fix some 7 > and consider the dynamics on the torus Tn generated 
by /^AT.pcr where the potential V satisfies Assumption [H Assume that the system is 
initially distributed according to a local Gibbs state / fTPj) with smooth energy profile Zq 
and volume profile Oq- Consider a positive time t such that the solution (e, o) to 
belongs to hi and is smooth on the time interval [0,t]. Then, for any continuous test 
functions G, if : T — )■ M, the following convergence in probability holds as N +00; 

{£N{tN,G)yN{tN,H)) ^ Q_^G{q)t{t,q)dq, j^H{q)^{t,q)dq 

The derivation of the hydro dynamic limits beyond the shocks for systems of 
conservation laws of dimension n > 2 is very difficult and is one of the most challenging 
problems in the field of hydrodynamic limits. The first difficulty is of course our poor 
understanding of the solutions to such systems. Recently, J. Fritz proposed in [13] to 
derive hydrodynamic limits for hyperbolic systems (in the case n = 2) by some extension 
of the compensated-compactness approach [2111101 [ID] to stochastic microscopic models. 
This program has been achieved in [15] (see also [E]), where the authors derive the 
classical n = 2 Leroux system of conservation laws. In fact, to be exact, only the 
convergence to the set of entropy solutions is proved, the question of uniqueness being 
left open. It remains nonetheless the best result available at this time. The proof is based 
on a strict control of entropy pairs at the microscopic level by the use of logarithmic 
Sobolev inequ ality estimates. It would be very interesting to extend these methods to 
systems such as the ones considered in this paper. 

4. Ergodicity 

We prove here the ergodicity of the stochastic dynamics, which is the fundamental 
ingredient for the hydrodynamic limit. 

4.I. Definitions and notation 

In order to explain what is meant by ergodicity of the infinite volume dynamics we 
need to introduce some notation. For any topological space X equipped with its Borel 
cr- algebra we denote by V{X) the convex set of probability measures on X. The relative 
entropy if(i^|yu) of z/ G V{X) with respect to /i G V{X) is defined as 

H{u\iJ.) = sup I [ (pdu -log ( [ d/i ) I , (27) 
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where the supremum is carried over all bounded measurable functions on X. Recall 
also the entropy inequality, which states that for every positive constant a > and every 
bounded measurable function 6, it holds 



(2^ 



Let 9xiX G Z, be the shift by x: {6x'r])z = rjx+z- For any function g on Q, OxQ 
is the function such that {6xg){ri) = g{6xri). For any probability measure fi G V{Q), 
6xfi G Vi^Q) is the probability measure such that, for any bounded function g : Q ^ M., 
it holds g d{9xf^) = 9xg dfi. If 9xf^ = for any x then /x is said to be translation 
invariant. 

If A is a finite subset of Z the marginal of G V{VL) on is denoted by /i|A- 
The relative entropy of G V{VL) with respect io ^ E V{VL) in the box A is defined by 
H{y\\ I /iIa) and is denoted by Hx{y\ix). We say that a translation invariant probability 
measure v G V{VL) has finite entropy density (with respect to /i) if there exists a finite 
positive constant C such that for any finite A C Z, H\{i'\^) < C|A|. In fact, if this 
condition is satisfied, then the limit 

^ lAKoo |A| 

exists and is finite (see [IS]). It is called the entropy density of u with respect to /i. 
We are now in position to define ergodicity. 

Definition 1. We say that the infinite volume dynamics with infinitesimal generator Q 
is ergodic if the following claim is true: If v E V{Q) is a probability measure invariant 
by translation, invariant by the dynamics generated by Q and with finite entropy density 
with respect to Hi^, then u is a mixture of the n^^x, /3 > 0, A G M. 

4-2. Ergodicity of the stochastic dynamics 

We are not able to prove the ergodicity of the deterministic dynamics in general, but we 
can prove it under the additional assumption that the invariant measure is exchangeable. 

Theorem 2. Assume that the potential V satisfies [24\)- f2^} . Let v be a translation 
invariant measure with a finite local entropy density w.r.t. fii^ such that 

V/GCo^(fi), jAfdu = 0. (29) 

If u is exchangeable then u is a mixture of fijs^x, /3 > 0, A G M. 

The proof of Theorem [2] is provided in Section 14.31 This result has an interesting 
consequence. 

Corollary 1. The infinite volume dynamics generated by C is ergodic. 

The proof of this corollary is similar to the proof given in [16] (or [3]). It is based on 
the fact that, if u is invariant for C, then it can be shown, by some entropy arguments, 
that u is invariant separately for A and for S. The invariance with respect to S implies 
that u is exchangeable and we can then apply Theorem |2] to conclude. 
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4-3. Proof of Theorem\^ 

We call J-inv the cr-field generated by the 6'i-invariant sets and J-^^^ the cr-field generated 
by the ^2-invariant sets. We denote by v the conditional measure i'{-\J^i 



;2 ) 

inv/ 



By the entropy inequality (!28|) . it is easy to show that 

Thus, the ergodic theorem gives the existence of the F^^^-vde&smakAe functions 

V^(r/) = lim ^ r/2x+i = 

\x\<e 



(30) 



a 

\x\<t 

where the convergence occurs in L^iy) and //-almost surely. The random variable 
can be considered as a constant under v. Let us fix x G Z and define 

/(r/) = (r/2.+i-V^)0(r/), 

where G C\{^ is a function depending only on the even sites {ri2z}zez- By Lemma El 
the conditional probability measure u is A invariant. Therefore, 

Afdu = 0= I (y'(r/2.+2) - V'{r]2.))(i>dv{r^) 

I (^'(^2.+i) - l^'(r/2.-i))(r/2.+i - V^)d^,.Adv 

+ j {V2.+1 - V) [V\v2.+3) - V'{v2.+i)] d,,,^^,,<l>du 

+ Yl jiV2.+l-V')[V'{v2z+l)-V'iv2z-l)]drj,Jdu. 

By Lemma [U, u is exchangeable. By exchanging 7722+1 and 7722-1 in the terms appearing 
in the last sum above, we see that the sum over z 7^ x, x + 1 is actually equal to 0. 
Moreover, by Lemma [1] again, the functions depending on even sites are independent of 
functions depending on odd sites under u, so that 



0= I {V'{r]2x+2)-V\r],,))^du 

(rr ^ ' (31) 

iV2x+i-V^)V'{r]2x-i)diy] / dn^^cpdu 



+ (1 iV2x+i-V')V'iv2x+3)du^ (1 d,,^^^,^<pdu^ . 
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By exchangeability of it holds, for any k ^ —1: 

y - v')v'{v2.-i)\K.) = y ((^2.+2.+i - v^)r(r/2.-i)| • 

Thus, 



(''72x+2fc+l 

|fc|<£,fc^-l 



7' 



The Ij^iv) limit in (130!) then gives 

i^((r?2x+i-Vi)y'(r/2._i)) =0. 
Similarly, it can be shown that z/((?72x+i ~ (?72x+3)) = 0. Since 

a^ = v [(r/i - V^) = v [(r/2.+i - V^) y'(r/2.+i)] , 

fl3T|) implies 

= 



(^'(^2(x+l)) - 1^'(r/2x))</'rfl^(r/) + " y (^'72. - •^^at.+i) j '^^i^' 

= y (r/2x+i - V^) dv. 



(32) 



In the same way, it can be shown that 

(^'(^2(x+i)+i) - y\ri2x^\))4>dv{r]) + a° 



(33) 



where = (j){{ri2z+i) zez) G C'o(^) is now a test function depending on the odd sites 
only, and V°, a° are the J-^^v measurable functions defined by 



a = u 



In fact, the 6'i-invariance of u gives 

= = z/(r/o| J-inv) := V, «° = «i = z/((r/o - V)V^'(r/o)| -Finv) := «• 

By Lemma [31 it holds a > //-almost surely. Taking into account the fact that {'r]2z)zez 
and {ri2z+i)z£Z are independent under z/, ( l32i) and ( l33i) allow to show, by Lemma HJ that 
the probability measure z/ is a product measure with marginals given by 

di^iVx ^ [r,r + dr]) = Z{V, a)^'^ exp (—aV{r) — Xr) dr, 

where Z[V, a) is a normalizing constant, and A := A(V, a) is such that 

Z(y, a)^^ / r exp {—aV{r) — Xr) dr = V. 



Let us summarize our results: Denoting by P the law of the random variables 
(a(cj), V(a;)) G (0, +oo) x M under i/, we have proved that, for any bounded local 
function /, 



This concludes the proof of Theorem |2l 
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Lemma 1. Let v he an exchangeable translation invariant measure with a finite entropy 
density w.r.t. fii Q. Then, z/(-) = is exchangeable and under v, the variables 

{ri2x)x& and {r]2x+i)x& are independent. 

Proof. For a given x G Z, consider the function T^'^+^ : Q ^ Q defined as (T^'^"*"^?]) = 
^x,x+i^ and denote by T be the set of local transformation T : Q ^ Q which are obtained 
as compositions of transformations T^'y~^^^ ?/ G Z. To show that u is exchangeable, we 
prove that, for a given x G Z and for any bounded function g, 

J go diJ = j g du. 

This amounts to proving that, for any J-^^^- measurable bounded function p, 

v{gp) = v(v {g\ p)=v{v{go T^^^+^\ j;^) p) = ^ [g o T^^^^'p) . 

By exchangeability of v., it is therefore sufficient to show that any J-^^^-measurable 
bounded function p is invariant by T^'^^^. 

To prove this result, we write p as an ergodic limit of local functions. Observe first 
that p = limfc^oo ^(Pfcl-^nv); ^ ^-S- and in L^ii'), where the local functions pk are defined 
as Pk = iy{p I J^aJ, with J^Afc the cx-algebra generated by {rj^ ; x G A^}. Besides, 

1 ^ 



so that 



1 ^ 

p = lim lim — — J2 ^2jPk- (34) 



Now, since pk is a local function, {62jPk) oT^'^+^ = 92jPk for \ j\ large enough. Therefore, 



-^x,x+l 



el e 



2£+ _ 



which, together with (IMl) . gives indeed p o _ 

We turn now to the second part of the lemma. Let 0° (resp. (p^) be a bounded 
local measurable function depending only on the even (resp. odd) sites. We have to 
show that 



^L) = H<t^'\^L) H<t^'\^L) ■ (35) 

Fix k,j G Z and consider a local transformation T & T such that cfp oT = 62k<P^ and 
(p^ oT = 92j(f)^. This is possible because 0° (resp. cj)^) depends only on the even (resp. 
odd) sites. Now, for any J^^^-measurable positive bounded function p, 

^V'l-^L) P) = ^ (0VV) = ((0° ° T) (01 o T) (p O T)) = (^2fc0° ^2,0' p) , 



1/ z/ 
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where the second equahty is a consequence of the exchangeabihty of i^, and the third 
one is obtained thanks to the invariance of p by T proved above. It follows that 



Taking first the limit £ — )■ oo and then the limit m — > cxo, we obtain, by the ergodic 
theorem, 



which is indeed (l35l) . □ 

Lemma 2. Let v a probability measure invariant by translation, invariant by A and 
with finite entropy density. Let f G Cq{Q) and p be a J-^^^ bounded measurable function. 
Then 

I {Af)pdu = u(^u{Af\TL)p) = 0. (36) 
In other words, i^(-|-7-^nv) is A-invariant. 

Proof. The tail cx-field J-'taii is defined by J^taii = ^k>iJ^Ai where J^\c is the a-field 
generated by {77^; |a;| > k}. Any J-^^^-measurable function p coincides u almost surely 



with a J-taii-nieasurable function (this can be seen for example as a consequence of ( 134|) ). 
In particular, for any /c G Z, ^^^.p = z/ a.s. . Thus, if / G Cq{Q), then pf G Co(f^) and 
by ( !29l) applied to the latter, we obtain = / A{pf) du = f (Af) p dv. □ 

Lemma 3. Let a = u ((770 — V)y(?7o)|Jinv) where V = z^(?7o|-^mv)- Then, 

u{{a < 0}) = 0. (37) 

Proof. Let u = i^(-|J-'inv) and remark first that a = i> ((770 — ?7i)l^'(?7o)), so that the 
exchangeability of u implies a = z/ ((771 — ?7o)l^'(?7i)), and 

a = ^HiVo-Vi)iV'ivo)-V'{vi)))- 

The convexity of V already gives a > 0. It remains to show that ^{{a = 0}) = 0. To 
prove this, it is sufficient to show that the restriction of i/|A to the box A = {0, 1}, has 
(z/ a.s.) a density with respect to the Lebesgue measure. In fact, we claim that the 
relative entropy of z>|a w.r.t. Po,i|a is finite, which implies the existence of the desired 
density. Indeed, consider any non-negative J-inv-measurable function p and any positive 
bounded J-A-measurable function ip. By the same argument as in the proof of Lemma [21 
p is J-taii- measurable. Consequently, by using this and the definition of the conditional 
expectation, we have 
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By the variational definition f l27p of tlie relative entropy and the previous equality, we 

get 



h{u\p, /iciU) + log e^rf/io,!^ ^. 



By assumption, v has a finite entropy density so that H 
the previous inequality is valid for any J-inv-measurable function p > 0, the claim 



Po,i|a) < +00. Since 



H V 



Po,i|a) < +00 follows by the variational formula of the relative entropy ( 127|) . 



□ 

Lemma 4. Let : M — )■ M 5e a 'potential satisfying Assumption^ and P > a constant. 
Let fi be a probability measure on Q satisfying, for any j G Z and any (p G Cq{Q), 

Then, fi is a product probability measure whose marginals are given by 
[Z{(3, X)]~^ exp{—l3V{rij) — Xrjj) drjj where A is such that v{f3,X) defined in /[T5\) coin- 
cides with fJ'irjj). 

Proof. Define ijje{ri) = e^^^i^^^^''*)"'"''*'''), where the value of A is chosen so that 

-l3Virjo)-Xvo 



-(svivo)-x-no 



Choosing 0(77) = x{v)i^i{v) with xiv) ^ local compactly supported smooth function, we 
obtain, for any j = !,...,£— 1, 

We now consider 

Xiv) = Xb{ri)g i^Yl ^» j ^o(^i' • • • ' 



where Xb is a local function not depending on rji, ... ,rie, and g is a. smooth function on 
M. Since 

if j = — 1, we can further condition on Ylk=iVk = and on the exterior 

configuration {rji, i 1, . . . ,i}, and obtain, for all j = 1, — 1, 



^Vk= ^u,7]i, i ^ 1, 



0.(38) 



k=l 
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These relations allow to show that the Borel measure 

ipe{v)f^ I dr]i, ...,dr]f: "^rjk = iu,r]i, !,...,£ 

\ k=l 

which has support on the hyperplane He = {{rji, . . . ,r]i) : Ylk=i = ^u}, is invariant by 
any translation of Tie. By Theorem 2.20 in [36], this measure is therefore the Lebesgue 
measure on Tie up to a multiplicative constant. This multiplicative constant is fixed by 
ensuring that the correct value of /i(?7o) is recovered, so that finally 

fi{dr]i, . . . , dr]i\r]i, i ^ I, . . . J) = dr]i . . . drji. 



Z{f],XY 



This concludes the proof of the lemma. 



□ 



4.4- Derivation of a triplet of compressible Euler equations from a perturbed Hamilton 
dynamics 

It is also natural to consider the Hamiltonian dynamics ([1]) perturbed by the random 
exchange of neighboring velocities p^ and Px+i, so that YlixP^^ Yli''^x and ^^.(t^(Px) + 
y{''^x)) are preserved during the time evolution. This has been considered in [16] in the 
case of a quadratic kinetic energy U{p) = p^/2. In the case of a non-quadratic kinetic 
energy U, it is possible to adapt the arguments of [16] to show that the corresponding 
dynamics is ergodic and that the triplet of compressible Euler equations ([5]) is obtained 
as a hydrodynamic limit by Yau's method. We give only a sketch of the proof of the 
ergodicity of the perturbed dynamics. 

Let U and V be two convex potentials satisfying Assumptions [H We state the 
dependence of the partition function with respect to the potential V (resp. U) by the 
notation Zy (resp. Zu). 

We consider the set f2 = f2 x equipped with its natural product topology and 
its Borel a- field. A typical configuration w G is denoted hj u = {r,p). Let 
{uj(t)}t>o = {r{t),p{t)}t>o be the the dynamics generated hj C = A + 71S, 7 > 0, 
with 

and S the momenta exchange noise operator, acting on test functions / : f2 x — )■ M 

as 

('5/)H = $^[/(r-,p-^'+^)-/(r,p)], 

x&Z 

where is the configuration obtained from p by exchanging the momentum p^ with 

Px+i- Observe that £, A and S conserve the energy, the momentum and the deformation. 
All the equilibrium Gibbs measures {1^/3, a,a' , /? > 0, A, A' G M} are invariant. 

Theorem 3. Assume that the potentials U, V satisfy Assumptions [3 The infinite 
volume dynamics generated by C is ergodic. 
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Proof. Let u he a. translation invariant probability measure with finite entropy density 
invariant for C We have to show this is a mixture of the z^/3,a,a'- As stated before, it 
can be shown by entropy arguments that u is separately invariant for A and for S. The 
invariance w.r.t. S implies that the law of p is exchangeable under u. Thus, we are 
reduced to proving that if u is invariant for A, shift invariant with finite entropy density 
and such that the law of p is exchangeable, then z/ is a mixture of the z/^,a,a''s- 

Let u as above and call J-inv the a- field of invariant events for the shift. As in 
Lemma [U or in the proof of Theorem 2.1 of [16], one can show that p and r are 
independent under = v{ 

' I inv ) • 

Moreover, since an invariant set is, up to a i^- negligible set, a tail invariant set, as 
argued in Lemma |2l it can be proved that v is invariant for A. Observe also that the 
law of p under v is exchangeable. Now, we can identify the distribution of r as follows. 
The thermodynamic relations give a one to one correspondence between the averages of 
the momentum and kinetic energy, (vr, T), and the chemical potentials, /3 > 0, A G M, 
through the relations 

7r(/3,A) = -9AlogZ^(/3,A), T(/3, A) = -9^ log A). 

Let /3, A be the J-inv-measurable functions defined by the relations 

7r(/3, A) = u{vk). T{f3, A) = u{U{pk)). (39) 

Consider a test function of the form ilj{r,p) = ip{r){pj — n) for a smooth function 
if. Since / Ail' dv = 0, by using the independence of r and p under u and the 
exchangeability of the law of p under z/, we obtain easily that 

J {V'{r,)-V\r,.,))ipdu + r' J {dr^^,^-dr^ip)du = (40) 

where /3~^ = J U'{pj){pj — 7T)du. The fact that (3~^ is well defined, i.e. /3 is stricly 
positive u a.s., can be proved as in Lemma [3l 

By Lemma HJ the equation (HOj) is sufficient to identify the distribution of r under 
u as a product measure Ylx&z^vii^^ ■^') 6xp{— /3l^(ra,) — AV^-jdra;. Moreover, f3 and A' 
satisfy 

Setting u = uirk) and using a test function of the form ip{p){rj — 
u), a similar computation shows that the law of p under u is in the form 
Ilx&z Zu0, A) exp{-f3U{p^) - \p^} dp^ where = / V'{rj){rj - u)du and 

-dxlog Zu0,^) = HPk), -d^log Zu0,X) = v{U{pk)). 

As above, /3 is z/ a.s. strictly positive. In view of (l39l) . we have (3 = [3 and A = A. 
Injecting these informations in the definition of we obtain /3 = /3. 

We have therefore proved that = J-'g xy and the conclusion easily follows. □ 
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5. (Sub) Diffusive scaling 

The hydrodynamic limit is nothing but a law of large numbers. The second step of the 
study consists in looking at the fluctuations. As explained in the introduction, we expect 
that fluctuations appear at a shorter time-scale than the diffusive one. The study of 
the fluctuation field ([9]) for the deterministic system is very difficult. A more tractable 
quantity which allows to decide whether a = 1 or a < 1 is the diffusivity 'Dp x \i{t). 

To define this quantity, we introduce some notation. The equilibrium 
compressibility matrix x is the symmetric matrix defined by: 

X= /i/3,A,A'(Po;i^o) /U/3,A,A'(Po;Po) /i/3,A,A' (Po; '^o) • (41) 

\ At/3,A,A'('^o;^^o) ^-p.xx'iro.po) fii3,x,x'{ro;ro) J 

Here, /i^,A,A'(/; fl') denotes the covariance of the two functions / and g w.r.t. fJ'i3,x,x'- 
Let Ix.x+i be the normalized current associated to the three conservation laws. We do 
not give a precise definition of the latter for the deterministic dynamics but the reader 
will translate easily in this case the one we give in (H7!) below for the energy-volume 
conserving dynamics. Roughly speaking, Ix,x+i is obtained from the usual microscopic 
currents of the three conserved quantities by a change of frame prescribed by the 
linearized flow (l8j) of the hydrodynamic equations. Let Cp^\^\i be the current-current 
correlation function defined as 

C^/3,A,A'W = $^( W(t)/0,l(0)* 

Then, the diffusivity is given by 



V^,x,x'{t) = x-' (l-^) Cp,x,x'{s)ds. (42) 

When the current- current correlation function is integrable, the limit when t — > +00 
of the diffusivity is well defined, and the following Green-Kubo formula is obtained: 

V'^^.^y = hm Vp,x,x'it) = / C^,x,x'is) ds. (43) 

The existence of the above limit depends on the time decay of Cp^x^y. For a diffusive 
(resp. super-diffusive) behavior, T>ji^\{t) is of order 0{1) (resp. of order 0{i}~")). 

The super-diffusive behavior of the deterministic dynamics (7 = 0) can be proved 
easily for the linear dynamics, i.e. V{r) = r^/2, and the Kac-van-Moerbecke dynamics, 
i.e. V{r) = + r — 1. In these cases, the dynamics is actually ballistic since the 
diffusivity is of order 0{t). For the linear dynamics, this follows from the fact that 
the total current is a constant of motion. For the Kac-van-Moerbecke dynamics, an 
application of the Mazur inequality (see |30]), as done in [35] and |3] for the Toda 
lattice, shows that the diffusivity is also of order 0{t). Apart from these two cases, 
showing a super- diffusive behavior for general potentials remains challenging. 

We now turn to the energy-volume conserving dynamics. In this case, we also do 
not know, in general, what the large time behavior of the diffusivity is. Nevertheless, 
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when l^(r) = we can compute explicitly the value of the current-current correlation 
function and deduce the behavior of the diffusivity. However, before stating precisely 
the result, we first need to modify the definition (H2l) of the diffusivity to account for 
the stochastic perturbation. We have now only two conserved quantities. The local 
energy-volume conservation is expressed by the formulas 



where the instantaneous currents f^Z+i^ fxl+i 
while the local martingales M|''^(t), M^''^(t) read 



(44) 



[\/(r/,+i(s-)) - F(r/.(s^))] d [iV,_i,,.(s) - 7^ 
[rix+i{s~) - ri^{s~)] d [A^^._i,^.(s) - 7s] , 



where (A^a;-i,a;(^))x6Z are independent Poisson processes of intensity 7. Denoting by 
[X, Y]t the quadratic variation at time t between two adapted processes X and Y , it is 
easy to see that the martingales MJ'^, M^'^ are such that 

[M^y\Ml% = Q, a,6e{e,t;}, y ^ z, 



m\M^^% = j / [Vivx+iis))-Vivxis))fds, 



[M:^\ M:^\ = 7 r [Vx+M - Ti,{s)f ds, 
Jo 

[M:'\M:% = ^ f [V{vx+i{s)) ~V{vx{s))][Vx+i{s) -Vx{s)]ds. 



(45) 



Let us also introduce the equilibrium compressibility matrix x'- 

-d/se -dpv 
-dxe -dxv 



X 



/i/3,A(V"(r/o); V"(r/o)) /i/3,A(^(r/o); ^0) 



(46) 



Observe that x is symmetric since d/^v = dxe. We define the normalized current Jx.x+i 



as 



Jx,x+l — Jx,x+1 ~ -ij(0 ~ {^x " ^ 5 Jx^x+l 



jx,x+l 



3: 



(47) 



x,x+l 



where the term -D^(^) (^j, — ^ has been subtracted in order to study fluctuations in the 
transport frame. We introduce now 



x,x+l 



it) 



Jx,x+l{s) - 7V 



ViVxis) 
Vx{s) 



ds + M2{t), M2{t) 
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Then, the diffusivity is defined by 

Let us rewrite this expression in a more convenient way, introducing he current- 
current correlation function 

Cp,x{t) = E ( Ai(0)* ) ^ . (49) 

Remark first that the terms in fHHj) coming from the discrete gradients in 'Wx,x+i 
disappear. Indeed, after summation over |x| < A^, only two boundary terms at x = ±A^ 
remain and since they are divided by 2N + 1, their contribution vanishes in the limit 
N — 7- +00. Besides, the cross terms between the martingales and the normalized currents 
Jy,y+i cancel. The argument is based on a time-reversal property of the current and can 
be read in [2]. Moreover, by P5|) . 



\ \\x\<N J \\x\<N J 



2t (2iV + l)7x. 



In conclusion, the diffusivity for the stochastic energy-volume conserving dynamics is 
given by 

T^lxit) = {^-'if CpAs)ds + 7 Id2, (50) 

where Id2 is the 2x2 identity matrix. Observe that the long time behavior of 'Dp^x{t) 
is clearly driven by the long time behavior of the current-current correlation function. 

Theorem 4. Define the function g : [0, 7r/2] x (0, +oo) M by 

g{uj,t) = le-8^*-"'(-) cos^u), 

and consider the infinite volume dynamics generated by C, with the potential V{r) = 
r^/2, started at equilibrium under the Gibbs measure /i/3,A- Then, 



y^(jx,x+i{t)Jo.iW , 

^-^ \ ' / /3,A 

x6Z 



/ g{u,t)duj 



(51) 



V 0/ 

It follows that the only non-trivial term of this matrix is of order as t goes to 

infinity. In particular T)'^ ^(t) is of order 0{t^^'^) . More precisely, 

Sy/t / (3 

""^'^^'^ ^ 3PV^ [ X 

Proof. To simplify the notation we omit the indices (3, A. We denote by ^ -, ■ ^ the 
semi-inner product defined on local integrable functions by 

«/,^?»=EK/(W)-(/)(^^)]=^lim^^^ J2 {0xf,9yg)-{f){g), 

x&Z °° \x\,\y\<N 
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and by "H the Hilbert space obtained by the completion of the vector space of local 
functions with respect to ^ ■, ■ Note that every discrete gradient, i.e. a local 
function of the form 9if — f, is equal to in "H. 

We introduce the Laplace transform matrix Ij{z) (for 2 > 0) of the current-current 
correlation function: 

POO 

Uz) = / e-'-' < Jo,i(s) , Jo,i(0)* > ds =< {z - £)"^Jo,i , J^,i >, 
where we adopted the short notation [z — Jo,i foi' 



(52) 



Since V{r) = rV2, it holds Z{(3,\) = (/3/27r)-V2e^V2^^ ^ = ^ = and 
e = ^ (AV/3 + !)• It follows that 

/ 1 A2 A \ 



X 



2/32 /33 /32 

A 1 



-2 y y - ^0 

Since any discrete gradient is equal to zero in "H, the only non-zero entry of L(2) is the 
(1,1) component 

L^'\z) =< {z - C)-%^i , r/0,1 >, 

where r]^^y = {rj^ - T){r]y - r). 

The determination of the current-current correlation function therefore amounts to 
solving the resolvent equation {z — C)u = r/0,1. Consider the vector space V spanned by 
the orthogonal basis {rix,x+k}x£Z,k>i (which is also the space spanned by the family fix,y 
for X ^ y). This space is stable by C since 

(a) Arfx^x+k = {0i — ^)Vx,x+k-i — {(^i — ^)Vx-i,x+k is the difference of two discrete gradients, 
hence is equal to in 

(b) when k > 2, Srfx^x+k = (1 + Oi)T]x^x+k-i + (1 + 0_i)r]x^x+k+i - 4^Vx,x+k, while 

^Vx,x+l = Vx-l,x+l + Vx,x+2 ~ 2?7x,x+l- 

We may thus look for a solution of the form 

= Pk{x)r]x,x+k, (53) 

xez k>i 

with the condition 

5^5^|p,(x)r<+oo (54) 

xez k>i 
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since u is sought in "H. In view of point (a) above, u should actually be a solution to 

oo 

x£Z k=l 

with Fi{x) = {z + 2'y)pi{x) — j{p2{x) + p2{x — 1)) and, for k > 2, 

Fk{x) = {z + A-i)pk{x) --i[pk-i{x) + pk-i{x + 1) + Pk+i{x) + Pk+i{x - 1) 

Using the fact that {rjx,x+k\x&,k>i is an orthogonal basis of V, and identifying the 
coefficients in front of the different terms, it follows that 



Fkix) 



-{k=l,x=0}- 



(55) 



Introducing the Fourier transform h{uj) (for a; G Ti) of a given function h € /^(Z, 



x6Z 



the conditions (l55l) can be equivalently reformulated as 

\z + 2^)pi{u) - 7(1 + e''^^)p2{uj) = 1, 

[z + 47)pfc(a;) - 7(1 + e-^'^^)pk.r{u) - 7(1 + e^'-^)pk+i{uj) = 0, k>2. 
By Parseval's relation, condition ( l54l) is equivalent to 



(56) 



/ |pfc(i^)|^'^t<; < +00. 



It is then easy to show that ( l56l) and the above integrability condition lead to 'pk(uj) 
Pi{uj){X{uj))''-\ with 



2 + 2;/(27) 



1 + e- 



2i7rti; 



1- A 1 



COS TTW 



1 + Z 



/(47)J 



and the boundary condition (z + 27)pi(a;) — 7(1 + e '^)pi{iu)X(uj) = 1. It follows that 



L''\z) =« u, 770,1 »= 5Zpi(a;) = r%{o) 



xez 



■ - 
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r 



where 



r(y) = [y + ^/a + yy-i . 

A simple computation shows that, for any y > 0, 



POO 


2 r/"^ 






'0 


TT io 



TT 



cos^fx) 



y + 2 sin (x' 



(ix = Till). 



Since the Laplace transform uniquely characterizes the underlying function, we deduce 
that the current-current correlation function is indeed given by (15T]) . 
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Now, by dominated convergence, the following limit holds as t — > +00: 
yt J g[u,t)du=— J exp I — 87tsm I — ^ M cos \ 



Jo V27r7 
We then obtain the desired result with ( 1501) . □ 

6. Steady-state nonequilibrium systems 
6.1. General setting 

The results of the previous section were limited to harmonic potentials (and the Kac- 
van-Moerbecke potential in the special case 7 = 0). For generic anharmonic potentials, 
we can only provide numerical evidence of the super-diffusivity. However, it is difficult 
to estimate numerically the time autocorrelation functions of the currents because of 
their expected long-time tails, and because statistical errors are very large (in relative 
value) when t is large. Also, for finite systems (the only ones we can simulate on a 
computer), the autocorrelation is generically exponentially decreasing for anharmonic 
potentials, and, to obtain meaningful results, the thermodynamic limit should be taken 
before the long-time limit. 

A more tenable approach consists in studying a nonequilibrium system in its steady- 
state. We consider a finite system of length 2A^ + 1 in contact with two thermostats 
which fix the value of the energy at the boundaries. The generator of the dynamics is 
given by 

'^7V,open = + l^N + XeB^jy^Ti + KBN,Try (57) 

where An (resp. Sn) is defined by ([ID]) (resp. ([H])) with An = {-N,...,N} and 
Bx,T = Td"^^ — V'{r]x)dr,^. The positive parameters Xi and are the intensities of the 
thermostats. 

The generator B^^t can be seen as a thermostatting mechanism since the semigroup 
{St)t>o generated by B^^t has, under suitable assumptions on V, a unique (reversible) 
invariant probability measure on M given by 

MdO = ZW, 0)e-/^^(«) rfe, /3 = T-i, (58) 

and Stf converges exponentially fast to ^'t(/) for any observable / G L^^iyr) (we will 
however not use these facts in the sequel). Observe that 

MV) = -dp{logZ{(3,0)) = e{f3,0). 

Hence, in order to fix the energy at site — (resp. A^) to the value ce (resp. e^), we 
have to choose I3e = T^^ (resp. = T~^) such that e(/3^,0) = (resp. e{l3r,^) = e^). 
In the special case when Ti = Tr, the Gibbs measure /^/j.oUiv is invariant, and it can 



then be shown (see Appendix A ) that the the law of the stochastic process associated 
to (!57|) converges exponentially fast to /i/3,o|Ajv- 
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The proof of the existence and uniqueness of an invariant measure in the case when 
Ti 7^ Tr is given in Proposition [1] for a suitable class of potentials (see Appendix A for 
a proof in the case 7 = and Appendix B for a proof when 7 > 0). We denote by (■)ss 
the unique stationary state for the dynamics generated by £Ar,opcn- 

Proposition 1. Assume that the smooth potential V satisfies 

• [Growth at infinity] there exist real constants k > 2, ak > 0, C > such that 

lim \~''V{\q) = ak\q\'^, 

, . . , (59) 

lim V'{Xq) = kak\q\ sign(g), 

A— >+oo 

V"(q) 

1™ = 0; (60) 

• [Non-degeneracy] For any g G M there exists m := m{q) > 2 such that 

V^'^\q)^Q. (61) 

Then, there exists a unique stationary probability measure for the Markov process 
generated by /^Ar^opcn- This stationary state has a smooth positive density with respect to 
the Lebesgue measure. 

The energy currents j^'^+i, which are such that CN,opcn(y{rix)) = —^jl^i,x (foi' 
X = —N, . . . ,N +1), are given by the expressions (jSj) for x = —N + 1, . . . , — 1 while 

£'?._i,_;v = A. [TeV'iv^j,] - {V'iv^j,))'] , 

J^3v+i = -K [TrV'ivN) - {V'ivN))'] . 

Since {CN,opcn(y{r]x)))ss = 0, it follows that, for any x = -N, ...,N + 1, (i^2+i)ss is 
equal to a constant Jj^{Ti,Tr) independent oi x. In fact, 

N-l 
x=-Af 

The latter equation is interesting from a numerical viewpoint since it allows to perform 
some spatial averaging, hence reducing the statistical error of the results. Finally, let 
us mention that, for finite systems, standard results of linear response theory (see 
for instance |4l [34]) allow to relate the nonequilibrium current J]^ to the current 
autocorrelation at equilibrium as 

where the expectation is over all initial conditions and realizations of the paths. 
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6.2. Harmonic potentials 

The current can be computed explicitly for harmonic potentials when no stochastic 
perturbation is present. In this ballistic behavior (no damping of the current 

with the system size) is found, as for standard harmonic oscillator chains [35]. It can 
also be shown that the macroscopic profile of energy g G [— 1, 1] — ?■ {V{ri[xq/N])) ss is flat 
and equal to (Ti + T,.)/4. 



Proposition 2. Consider the system described by the generator (dl) with 7 = 0, for 
the harmonic potential V{r) = Then the steady-state energy current is 

Proof. Introduce Cx,y = {'r]xVy)ss = Cy^x for (x, G A^. First, note that the stationarity 
of the current implies Cx,x+i = — -^Ar(T£, Tr) for x = —N, . . . ,N and 

^e(Ti — C_Ar_Ar) = —Xr(Tr — Cn,n) = <^Ar , ^r) • (64) 

Besides, for x = -N + 1, . . . , - 2, it holds 

('^7V,opon(^a;^a;+l))ss = = Cx+l^x+1 + C'x,x+2 ~ (C*x,x + Cx-l^x+l)- 

Therefore, C7v_i,Ar_i + Cn-2,n = C^n+i-n+i + C.n-n+2- Now, 

('CAr,open('7-A''7-A^+l))ss = = C^N+1~N+1 + C^N+2,N — XiC^N,-N+l — C-N-N, 

and, similarly, Cn,n-2 + Cn-i,n-i = Cn,n — KCn,n-i- In conclusion, 

-N-N+l + C -N-N — Cn,N — KCn,N-1- 



The expression (1631) is then obtained by combining the above equality with f l64l) and the 
relation Cn^n-i = C_n-n+i = -J%{Ti,Tr). □ 

When 7 > 0, we expect to observe an anomalous diffusion with divergence rate 
5 = 1/2, i.e. NJ]^ ~ C^^/N when N is large enough. This is indeed confirmed by 
numerical simulations, see Figure [TJ 

6.3. Anharmonic potentials 

The nonlinear case is much more difficult. We estimate the exponent S > such that 

NJjj ~ A^-^ (65) 

using numerical simulations. If 5 = 0, the system is a normal conductor of energy. If 
on the other hand 5 > 0, it is a superconductor. 



30 



6.3.1. Numerical scheme. The time-discretization of the dynamics with generator fl57l) 
is done with a standard sphtting strategy, decomposing the generator as the sum of 
a deterministic part An, a thermostat part \iB-n,Ti + KBN,Try and the stochastic 
perturbation 'ySjsf, and integrating each part in this order. We denote by At the time- 
step. 

A simple numerical scheme for the deterministic part of the dynamics by relying on 
the Hamiltonian interpretation of the system. The longtime integration of Hamiltonian 
system is well understood. The most standard scheme used in practice is the so-called 
Stormer-Verlet scheme (see again [20]), which, for separable Hamiltonians (where the 
total energy is the sum of a kinetic part depending only on the momenta, and a potential 
part depending only on the positions) can be seen as a Strang approximation of the 
Hamiltonian evolution with positions and momenta updated successively. This amounts 
here to updating successively variables with odd and even indices. This corresponds to 
the Strang splitting based on the following decomposition of the generator: 

/I /I even , /lodd 

JiN — -^N + -^N 5 

with, in the case when is even. 



^even ^ ^ (v'{v2.+1-n) ' V {^2.-1- n)) 0,,^^^ + V {v.N+l)d,_, ' V\m-l)d,^, 



and a similar definition for A^"^. This splitting is particularly convenient since the time 
evolutions generated by are A'}^'^ are both analytically integrable. In conclusion, 

the numerical scheme used for the deterministic part reads: 



: + (^'(C+i) - V\r^:.,)) , X = -iV + 1, -iV + 3, . . . , AT - 1 
^T" = Vy + At {v'iv'^T) - V'iv;T)) , y = -iV, -iV + 2, . . . , AT 

V:^' = r,:^''' + ^ {v\r,:Xl) - V\r,:tl)) , X = -AT + 1, . . . , AT - 1. 

The thermostat part is taken care of by a simple Euler-Maruyama discretization: 



'r]%^ = V-N - XiAtV'{'nl^) + ^/2XeTeAtG'_ 

= Vn- A.At V^'(7/]^) + v/2A,r,AtG^, 

where {G'±^) are independent and identically distributed standard Gaussian random 
variables. Finally, the noise term with generator 71S is simulated by exchanging 773. 
and rjx+i (for x = —N, . . . ,N — 1) at exponentially distributed random times, with an 
average time 7"^ between two such exchanges. More precisely, we attach to each couple 
(x, X -|- 1) a random time r^, with r° drawn from an exponential law with parameter 7, 
and where {Tl")i^m are independent. This time is updated as follows: If > At, then 
T^~^^ = r™ — At, otherwise rj^ and rj^+i are exchanged and a new exchange time t^~^^ 
is resampled from an exponential law of parameter 7. 
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Figure 1. Current Jj^ as a function of the system size 2N + 1 for harmonic potentials 
V{r) = (in logj — logj scale), with three values of 7: 7 = (black), 7 = 0.1 (blue) 
and 7 = 1 (red). 



6.3.2. Numerical results. We considered the following potentials: 

(i) harmonic potential V{r) = 

(ii) anharmonic FPU-like potential V{r) = + 

(iii) Kac-van-Moerbecke potential V{r) = e"'' + r — 1; 

(iv) rotor V{r) = 1 — cos(r). 

The time-step At is chosen to ensure a good longtime preservation of energy for the 
deterministic dynamics in the absence of stochastic perturbation. In order to have a 
relative error in energy less than 10^^, we used At = 0.005 except for rotors where 
At = 0.025. We set = = 1 (this choice maximizes the observed current in 
the harmonic case according to fl63l) ). and considered a small temperature difference 
Ti = 1.1, Tj. = 0.9. We performed 10^ iterations for < 2^^, 5 x 10® iterations for 
= 2^^, 2.5 X 10® iterations for = 2^^ and 1.25 x 10® iterations for = 2^^ The 
corresponding system size in the latter case is 2A^ + 1 = 65, 537. 

The simulation results are presented in Figures [1] to HI The conductivity exponents 
extracted from the numerical simulations presented in Figures [1] to [3] are reported in 
Table [H Exponents in the harmonic case agree with their expected values (see the 
discussion in Section [6l2l) . For nonlinear potentials, except for the singular value 6 = 1 
when 7 = 0, the exponents seem to be monotonically increasing with 7. We conjecture 
that these exponents should attain the limiting value 0.5 as 7 — )■ +00. A similar behavior 
of the exponents is observed for Toda chains [22]. Note also that the value found for 
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10 11 12 13 
log_2(size) 



Figure 2. Current as a function of the system size 2N + 1 for anharmonic potentials 
V{r) = r^/2 + r*/4 (in logj — log2 scale), with four values of 7: 7 = (black), 7 = 0.01 
(blue), 7 = 0.1 (pink) and 7 = 1 (red). 
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Figure 3. Current as a function of the system size 2N + 1 for the Kac-van-Moerbecke 
potential V{r) = e^'' + r — 1 (in logj — logj scale), with five values of 7: 7 = (black), 
7 = 0.001 (blue), 7 = 0.01 (green), 7 = 0.1 (pink) and 7 = 1 (red). 
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Figure 4. Current as a function of the system size 2N + 1 for a chain of rotors: 
V{r) = 1 — cos(r). A normal conductivity is observed since the estimated conductivity 
exponent is 5 ~ 0.02. 



Table 1. Conductivity exponents 5 (see (|65|) for the definition of S). 

7 harmonic anharnionic KVM 

1 0.13 1 
0.01 - 0.14 0.12 
0.1 0.50 0.27 0.25 

1 0.50 0.43 0.33 



7 = with the anharmonic FPU potential V{r) = + is smaller than the 
corresponding value for standard oscillator chains, which is around 0.33 (see [2H])- 



Appendix A. Existence and uniqueness of the stationary state for 7 = 



We adapt in our context the methods introduced in [5] (see [M] for a review) for 
the chains of coupled oscillators, following closely the exposition given in [7]. The 
proof is divided in 3 steps: (i) the aim of Appendix A. 2 is to prove the Lyapunov 
condition f lA.SP : (ii) we then prove the smoothness of the transition probability using 



hypoellipticity arguments (Appendix A.3), and finally (iii) show that the process is 



irreducible ( Appendix A.4 ). These three arguments, together with [33] Theorem 8.9] 
allow to conclude the existence and the uniqueness of an invariant measure with smooth 
density, as well as the exponential convergence of the law of the process to the invariant 
measure. 
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Appendix A.l. Definition of the process 

The process generated by Cn, open is denoted by {rj{t))t>Q and is the solution of the system 

Ur^^Nit) = V'{r]^N+i)dt - XeV'{r]^N)dt + ^2\iT,dB^Nit), 

(S) \ dr]^{t) = (y'(r/..+i) - V'{v.-i))dt, x = -N + 1, . . . , N - 1, 

[dr]Nit) = -V'{7]N-i)dt - XrV'{r]N)dt + ^/2KT^. dB^it), 

where B^n, Bn are two standard independent Brownian motions. Since the drift 
coefficients are only locally Lipschitz, it is unclear whether a solution to (S) exists 
for any time. This problem can however easily be solved by using a suitable Lyapunov 
function, as we now show. For a > 0, consider 



N 



Wa{v)=exp\a V{v.)\- (A.l) 

x=-N 

A simple computation shows that 



+ XeTeV"iv-N) + MT,a - l)iy\ri.N)) 



so that, using the assumption (l60l) . CN,o-pc-nWa < AWa for some positive constant A 
as soon as a is sufficiently small. This is sufficient to prove the existence of a unique 
solution to (S) for any initial condition and for any time t > 0. We denote by {Tt)t>o 
the corresponding semigroup. 

Appendix A. 2. Proof of the Lyapunov condition M.5]) 

Given a solution {rj{t))t>Q of (S), we define the scaled process rj^ (for E > Q) as follows: 

ri^it) = E-^'^'riitE^/^-^). (A.2) 

The scaling is chosen so that there is no explicit dependence in E in the formal limit 
E — 7- +00 (see the system (Soo) below). The process rj^ is a solution of 



'dr^%{t) = E'/'-'{v'iE'/'v^^^,) - X,V'iE'/'v^j,))dt 



+ xl^dB.^{t), 



(Se) < 



dv^it) = i?^/^-^(v'(i?^/Sf+i) - V'iE'/'^r^l, 

x = -N + 1,...,N -1, 
dvm = -E'l'-'{y\E'IS''N_,) + XrV\E'/'r,f,))dt 



dt, 

(A.3) 



+ J^-^dB^it), 
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where -B-tv, are two standard independent Brownian motions. Observe that the 
scaled process is such that if i7(?7(0)) = E then HE{ri^{0)) = 1, where the scaled energy 
He is 



x=-N 



In the limit E — )• +oo, the noise disappears and, by the scaling property (!59l) of the 
potential, (Se) reduces formally to the following deterministic system: 



(Soo) < dr],{t) = (^e'ivx+i) - ^'(r/._i))rft, 

dr]N{t) = -(d'iVN-i) + \rO'{r]N)^dt, 



X 



where 6{q) = ak\q\'' is a function since k > 2. We also introduce the corresponding 
limiting energy function 

N 



-N 



Arguing by contradiction, it is easy to show that if {ri(t))t>o is a solution of (Sc 
starting from an initial condition rj{0) such that iJoo(?7(0)) = 1, then, for any r > 0, 



[{e'iv-NisW + {o'ivNisW] ds > 0. 



(A.4) 



By using the continuity of solutions of stochastic differential equations with respect to 
both parameters and starting points, we can then state the following asymptotic result. 

Lemma 5. Assume that {{r]^{t))t>o)nm is a sequence of solutions of (S) starting from 
?7"(0), with En = i/(?7"(0)) — t- +oo. Then there exists a subsequence ((77™(t))t>o)meN 
such that, for any C > and to > 0, 

-to 



lim E, 



exp 



-C 



[iO'iV-Nis)r + (O'iv'^^is))'] ds 



0. 



Proof. Recall that {ri^'^"{t))t>o is the solution of (S_b„) starting from 77'^'-^" (0), which is 
such that HE^{f]"'^''{0)) = 1. This implies that the sequence {f]"''^"{0))n£N remains in 
a compact set and we can extract a subsequence, denoted by (?7'^'^'"(0))meN, such that 



m,E„ 



(0) converges to some element rj* G 



-N,...,N} 



By continuity, Hooij]*) = 1- Let 



{rj{t))t>Q be the solution to (Sqo) starting from rj* . 

We fix C > and to > 0. For any r > 0, thanks to the continuity of solutions of 
stochastic differential equations with respect to both parameters and starting points, it 
holds 



lim E 

m— >oo 



sup |?7 

t<T 



rri,E„ 



0. 
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For any 77 G M''^ n,...,n}^ denote the term [(^'(?7_Ar))^ + (^'(?7Ar))^] by Kirf). Then, for 
any a > 0, 



hm E 

m—^oo 



exp (-a / i^(r/'^-^'"(s))rfs 



-a/(r) 



where /(r) = K {ri{s)) c/s > by (]A.4p . Observe also that, for > 1, 



to 



K{r]^{s))ds = E„ 



Ki^q'^'^^i^u)) du>Em K{r]"'^''^{u)) du. 



to 



Now, consider some arbitrary constant A > 0. By choosing m sufficiently large so that 



A < CEfn, it follows 



lim E, 



The result is then obtained by letting A go to infinity. 



exp (^-C K{7]'^{s))ds^ 



□ 



We now claim that there exist a time to > 0, finite constants 6„ and k„ G (0, 1) 
with lim^^oo i^n = 0, and compact sets Kn such that 

where := Wa is the Lyapunov function ( lA.ip . By choosing compact sets Kn of the 
form {t] G | W{r]) < an} with a„ — >■ +00, it is enough to show that 

hm sup = 0. 

W{7]) 

Therefore, (lA.Sp is a consequence of Lemma [5] and the following result. 

Lemma 6. // a is sufficiently small, there exist c, C > and q > 1 such that, for all 
r/(0) G Ri-^'-'^i, 



TtoWim) 



< e'*" E. 



''?(0) 



to 



exp -C / [(^'(r/_;v(s)))' + {0'iVNis)))'] ds 



w{m) 

Proof. It holds 

H{r^{t)) = H{r,m+ ! Cr,,openH{r]{s)) ds + M{t), 

Jo 

where M{t) is a continuous martingale of quadratic variation 

[M], =^ (/:^,opcn(i^')-2i/£^,openi^)(r/(s))rfs 

A,T,[V'(r/_jv(s))]=^ + XrTr[V'{Ms))f)ds. 
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Now, for any constants p,q > 1 such that p ^ + q -'^ = 1, 
TtW{r]{0)) _ r a{H{v{t))-H{v{om 

= E^(o) exp(aM{t) + a {CN,opcnH){r]{s))ds 







'»?{0) 



exp \aM{t) - p— [M], + p— [M], + aj {CN,oponH){r]{s))ds 



where 



E,(o) [XtYt] < E,(o)[Xf]i/''E,(o) [Yt' 



XI = exp ( paM{t) - [M], 



2 

is an exponential martingale with constant mean equal to 1. Moreover, 
Fj*^ = exp (^PQy l^]t + {CN,opcnH){r]{s))ds^ 

= exp (^j (^F{r]^N{s), Xe,Ti) + F{r]N{s), Xr,Tr)jds^ , 

with 

F(r, A, T) = Aag [(paT - l){V'{r)y + Ty"(r)] . 

Taking a sufficiently small, and using (159|) and (l60l) . we see that there exists two 
constants C, c > (depending on a,p, T^, T^, A^, A,, and V) such that 

V(r, A, T) G M X {A,, AJ x {T,, TJ F(r, A, T) < gc - C[^'(r)]2. 

This completes the proof of the lemma. □ 

Appendix A. 3. Smoothness of the transition probability 
The generator £Ar,open can be written as 

jOpen 

where Xq,X±i\j are ffist-order differential operators: 

Xo = An- XeV'{r]^N)dn_j^ - XrV'{riN)d^j,, X±Ar = 9^^^. 

If the Lie algebra £ generated by the vector fields Xq, X±n, i.e. the smallest Lie algebra 
containing 

has full rank at every point 77, then the generator £Ar,open is hypoelliptic, and, 
by Hormander's theorem on hypoelliptic operators (see [21]), the semigroup {Tt)t>Q 
generated by /^Ar^opcn has a smooth transition probability density, is strong Feller and 
the invariant measures, if they exist, also have a smooth density. 
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To show that £ has full rank, we first observe that G £. We then prove that 
d.^^^^ e £. Indeed, [Xo,r]N] = -V" {riN)d^^^^ + XrV"{riN)d^^, so that -V"{riN)d^^_^ e 
2. If V"{riN) 7^ we already have 9,,^ ^ G £. Otherwise we compute the iterated 
Lie bracket [. . . , [Xq, dr,j^], . . . , 9^^^] and obtain V^'^\r]N)drij^_^ G -C for any m > 2. The 
non-degeneracy condition ( I6T1) on V therefore gives drjj^^-^ G £. By iterating the above 
argument, it follows easily that G £ for every x G {— iV, . . . , A^}. 

Appendix A. 4- Irreducibility of the dynamics 

We show here that the semigroup {Tt)t>o is strongly irreducible, that is, for every 77, 
every t > 0, and every non-empty open set A, it holds Tt{ri,A) > 0. Recall that, in 
general, irreducibility is not a consequence of hypoellipticity. 

Irreducibility can be proved by using a well-known relationship between stochastic 
differential equations and control theory (see for instance the discussions in [9l Section 3] 
or [71 Section 4]). In [19], a general approach for obtaining the desired controllability 
is presented for divergence-free systems having a conserved quantity and satisfying 
a Hormander condition. Since our system satisfies these assumptions, we can apply 
Theorem 2.1 of [19] and deduce that the semigroup (Tt)t>Q is strongly irreducible. 



Appendix B. Existence and uniqueness of the stationary state for 7 > 

We denote by (Tj)t>o the semigroup generated by ^Ar^opcn for some given 7 > 0, and 
by (Tt)t>o the semigroup corresponding to CN,opcn for 7 = 0. The same arguments as 
for the case 7 = can be used to show that (Tf)t>o satisfies some Lyapunov condition 
similar to ( ]A.5I) . Relying on [331 Theorem 8.1] for instance, it is then enough to show 
that the process is irreducible and that the transition probability has a smooth density. 



Appendix B.l. Irreducibility of the dynamics 

We show first that {Tt)t>o is strongly irreducible. Let us denote the probability transition 
of Tt hj pti^r],^) and the probability transition of Tt by pt{rj,d^). 

Let 0"! be the stopping time defined as the first time when two variables rj^^ and 
r]x+i are exchanged. Observe that o"i has an exponential law of parameter 277V. For 
every bounded measurable function / : Mi^^' - '^i it holds 

(Tl<tJ 

h (B.l) 

+j fdse-'-'"^ ^ [d^Ps{v,o( [ pt-sr^''^\dafia 

Jo Ji \Ji' 



39 



Iterating the above argument, we obtain the following formula for pf. 



fc=l xi,...,Xk=-N '-"'0 



This shows that Pt{v,d^) = Pt{fliC)d^ is absolutely continuous with respect to the 
Lebesgue measure. Moreover, the semigroup {Tt)t>Q is strongly irreducible because 
Tt is strongly irreducible. 

Lemma 7. The semigroup {Tt)t>Q is strongly Feller, i.e. it maps bounded measurable 
functions to continuous bounded functions. 

Proof. The semigroup {Tt)t>o is strongly Feller. This implies (see e.g. [381 Corollary 2.4]) 
that for every t > and every compact set K, 



lim sup sup ( Ttu ) (v) - { Ttu ) {r]'] 

<5^0 \V-V'\<S ll«l|oo<l 
ri,rj'eK 



0. 



(B.2) 



Let / be a bounded measurable function with ||/||oo < 1- We have to show that, for 
any fixed t > 0, Ttf is a continuous bounded function. By (IB.ip we have 



N-l 



+ 7 E 



x=-N 



'2yN{t-s) 



Tt.s o o T, o / (r/) - T,_, o o T, o / (r/') ds. 



Here, is the operator acting on functions / = /(r/) as {Fxf){ri) = f{ri^'^^^). Observe 
that the absolute value of the second term on the right hand side is bounded above by 



x=-N 



llffl|tx><l 



Tt-sQ )iv)- ( Tt^sQ iv') 



ds. 



By the bounded convergence theorem and (IB. 21) we have 
lim(T,/)(r/')-(T,/)(r/) = 0, 

which concludes the proof. 



□ 



These two last properties (irreducibility and strong Feller property) are sufficient 
to have uniqueness of the invariant measure fiss- To show that the latter has a density, 
we observe that for any t > 0, the condition UssTt = Tt implies that, for any measurable 

set A of M{-^.-..,Af} 
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where the second hne follows from Fubini's theorem. 
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